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Research Article 
ABSTRACT 


The outbreak of COVID-19 affects millions of people across the world. 
Challenges and burdens were faced in the healthcare system due to the 
absence of effective antiviral drugs. Herein, we attempted the repurposing of 
the FDA approved antiviral drugs from both synthetic and marine sources 
against the SARS-CoV-2 Main protease (Mpro). Interrupting the enzymatic 
activity of Mpro is a feasible approach for drug development against the virus 
as it plays an integral role in the life cycle of the virus by its involvement in 
viral replication and infection. This study attempts to identify a potentially 
interacting drug molecule against Mpro by molecular interaction analysis and 
validation. The efficacy of the antiviral drug molecules from synthetic and 
marine sources was compared using multiple docking software. The interaction 
was found to be stable after MD simulation and MM-PBSA analysis. 
Pseudopterosin A and Simeprevir show a noticeable binding affinity of -7.9 
kcal/mol and -8 kcal/mol respectively. The stable interaction after MD 
simulation and better binding free energy indicates the significant potency of 
the molecules to bind the active site of SARS-Cov-2 Mpro. The study demands 
further clinical validation of the resultant drug molecule to act against the 
virus 
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INTRODUCTION 


SARS-CoV-2 is single stranded RNA enveloped viruses consisting of structural, accessory, and non-structural proteins. 
Envelope (E), Spike (S), Nucleocapsid (N), and Membrane (M) is the structural proteins whereas RNA dependent RNA 
polymerase (RdRp), 3-chymotrypsin-like/main protease- (3CLpro/Mpro), and Papain Like protease (PLpro) are the non- 
structural proteins. Mpro and PLpro have a fundamental role in viral replication as they synthesize viral polyproteins 
via multiple mechanisms [4]. 

Replication of the virus in the host cell is carried out by replicase polyprotein, which is facilitated by the Mpro. The 
Mpro has three domains; domain-1 (residues from 8-101), domain-2 (residues from 102-184), and domain-3 
(residues from 201-303) in which an antiparallel B-barrel structure confers the first two domains, and the third 
domain consists of five a-helices. Among all coronaviruses, Mpros share a structurally highly conserved substrate- 
recognition pocket and catalytically functions as a homodimer (protomer A and protomer B), each protomer comprising 
three domains. The functional significance and the lack of homologous proteins in humans make Mpro an attractive 
target for the drugs. Interrupting the catalytic activity of Mpro is a relevant strategy for anti-Covid drug development 
(2, 


Recently, drug repurposing plays a significant role in expediting drug development against exponentially expanding 
diseases. Repurposing drugs is an approach for figuring out a novel application for approved or clinical drugs beyond 
the scope of the primary medical purpose. The complexity and convolution of the normal drug development process 
and the necessity of drugs for the pandemic direct us to choose the repurposing of existing drugs. Remdesivir, 
hydroxychloroquine, lopinavir-ritonavir, ivermectin, etc, were repurposed for the medication against SARS-CoV-2. 

The unique oceanic habitat provided a wide variety of natural resources with important therapeutic activities like 
antiviral, antioxidant, antimicrobial, etc. Vidarabine, an extensively used antiviral drug for the therapy of HSV 
infection, was derived from spongouridine, a secondary metabolite from a marine sponge. Portimine and rhaman 
sulfates, derivatives from Vulcanodinium rugosum and Monostroma nitidum respectively, had reported anti-viral 
activity against HIV and influenza virus. The mentioned viruses are single stranded RNA viruses like SARS CoV-2 hence 
the marine based drugs can be exploited for the medication of the same !°1. 

The SARS-CoV-2 virus has similar phylogenetic features to the earlier versions; MERS-CoV and SARS-CoV, hence 
therapy using antivirals can be considered as a guide to developing novel therapeutics for the treatment of COVID-19. 
The marine microbes are a proven source of anti-infective agents from which various drug lead molecules were 
extracted. FDA approved marine-based drugs with antiviral properties and other synthetic proven antiviral drugs were 
considered for the study. Here, we have repurposed these drug molecules to find their interaction efficacy against the 
Mpro of the SARS-CoV-2 virus (4). 


MATERIALS AND METHODS 


Target structure selection and retrieval 

The 3D structure of SARS-CoV-2 main protease (PDB ID: 6LU7, the crystal structure of COVID-19 main protease in 
complex with inhibitor N3) with a resolution of 2.16 A was retrieved from the protein data bank. The physical-chemical 
properties and secondary structure prediction of the selected protein structure were computed using Protparam and 
SOPMA respectively. Single chain-A with 306 amino acid residues in the target structure was considered for the study. 
The water molecules, unwanted chains, and attached ligands were removed using the Discovery Studio visualizer 
followed by adding missing residue using the Swiss Pdb viewer software !51. 
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Table 1. Physicochemical properties of 6LU7. 


PDB /T |Al | MW No. ofamino | GV | Il Half-life Extinction coefficient 
ID Pl acid Mammal Cys 
ian Yeast | E.coli | Cys reduced 
5.9 | 82.1 | 33796. 0.01 | 27.6 >20 >10 33640/0. | 32890/ 
6LU7 | 5 2 64 306 9 5 1.9 hrs hrs hrs 995 0.973 
*T-Pl: Theoretical Pl; “Al: Aliphatic Index; “MW: Molecular Weight; “GV: Gravy Value; “Il: Instability Index 


Table 2. Secondary structure prediction of 6LU7. 
Alpha helix percent 
29.08% 


PDB _ID 
6LU7 


Extended strand 
27.12% 


Random coil 
32.35% 


Beta turn percent 
11.44% 


Ligands 

A set of 40 commercially available antiviral drugs (marine-20, synthetic-20) were considered as ligand molecules for 
the study (Table 3). The molecules were screened and listed for their anti-viral activity through the literature survey. 
And their 2D structures were retrieved from the PubChem database in sdf format. The structures were converted to 
pdbqt format using the Open babel tool !6!. 


Table 3. The ligand molecules selected for the study (20 marine and 20 synthetic) with their source, application, and 


CAS number. 
S.No Drug name Source Application CAS no 
Used to treat herpes simplex 
Marine (Cryptotethya virus infections, chickenpox, and | 59277-89- 
1 Acyclovir crypta) shingles. 3 
HIV-1 synthesis of the natural 
UAG suppressor glutamine 
transfer trna; HIV-1 crossing the | 55303-98- 
2 Avarol Marine (Disidea avara) blood brain barrier 5 
Marine (sponge, Tethya 30516-87- 
3 Azidothymidine cripta) Used to treat HSV, HIV infection 1 
Used to treat susceptible 
infections with P. Vivax, P. 
Malariae, P. Ovale, and P. 
4 Chloroquine Marine (C. vulpina) Falciparum. 54-05-7 
Used to treat acute non- 
lymphocytic leukemia, 
lymphocytic leukemia, and the 
Maine (Cryptotethya blast phase of chronic 
5 Cytarabine crypta) myelocytic leukemia. 69-74-9 
Marine (Luffariel! 75088-80- 
6 Manoalide variabilis) Anti-inflammatory sesterterpene 1 
HIV-1 RNA and DNA dependent 
DNA polymerase activities of 
7 Polyacetylenetriol Marine (Petrosia sp.) retro viral RT - 
Marine 
(Pseudopterogoria 104855- 
8 Pseudopterosin A Elisabethae) Antiviral drug 20-1 
Marine 
(Pseudopterogoria 104855- 
9 Pseudopterosin B Elisabethae) Antiviral drug 21-2 
Marine 
(Pseudopterogoria 104881- 
10 Pseudopterosin C Elisabethae) Antiviral Drug 78-9 
Marine 
(Pseudopterogoria 104855- 
11 Pseudopterosin D Elisabethae) Antiviral drug 22-3 
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Marine 
(Pseudopterogoria 121011- 
12 Pseudopterosin E Elisabethae) Antiviral drug 80-1 
Marine 
(Pseudopterogoria 
13 Pseudopterosin F Elisabethae) Antiviral drug - 
Marine 
(Pseudopterogoria 
14 Pseudopterosin G Elisabethae) Antiviral drug - 
Marine 
(Pseudopterogoria 
15 PseudopterosinW Elisabethae) Antiviral drug - 
Marine 
(Pseudopterogoria 
16 Pseudopterosin X Elisabethae) Antiviral drug - 
Sarcophine derivatives are used 
Marine (Sarcophyton as anticancer and anti- 55038-27- 
17 Sarcophine glaucum) inflammatory agents 2 
72629-69- 
Marine (Sarcophyton T/ 
18 Sarcophytol glaucum) Cancer chemopreventive agent 5284454 
Sarcophine derivatives are used 
Marine (Sarcophyton as anticancer and anti- 
19 Sarcophytolide glaucum) inflammatory agents 11381519 
Antiviral medicine used to treat 24356-66- 
20 Vidarabine Marine (Tethya crypta) the herpes simplex virus. 9/ 
Synthetic (nucleotide Used to treat (chronic) infections 106941- 
24 Adefovir analogs) with hepatitis B virus 25-7 
Synthetic (synthetic 
derivative of 
hydroxyethylamine 161814- 
22 Amprenavir sulfonamide) Used to treat HIV infection 49-9 
Synthetic (quinazoline 68475-42- 
23 Anagrelide derivative) Thrombocythemia 3 
Synthetic (nucleoside 
analogue of 2’- 142217- 
24 Entecavir deoxyguanosine) Treatment of hepatitis B virus 69-4 
Synthetic (purine Treatment of viral infections 259793- 
25 Favipiravir nucleic acid analog) including influenza. 96-9 
Synthetic 
(pyrophosphate Treat herpes simplex virus, to 
26 Foscarnet analogue) treat cytomegalovirus 4428-95-99 
Hepatitis C NS3/4A protease 
Synthetic (second- inhibitor used to treat hepatitis 1365970- 
27 Glecaprevir generation Pl) C. 03-1 
Synthetic (derivative of Used to prevent and treat 
28 Hydroxychloroquine chloroquine) malaria 118-42-3 
Synthetic (N-(2- 
hydroxyethyl)piperazine, 
a 
piperazinecarboxamide, 
and a dicarboxylic acid Treatment of patients with HIV 150378- 
29 Indinavir diamide) infection (Pl) 17-9 
Synthetic 
(allophenylnorstatine- 186538- 
30 JE-2147 containing dipeptide) HIV protease inhibitor 00-1 
Synthetic (analogous 216863- 
31 L-756423 with Indinavir) HIV protease inhibitor 66-0 
32 Lopinavir Synthetic (dicarboxylic Treatment of patients with HIV 192725- 
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acid diamide) infection (Pl) 17-0 
Synthetic 
33 Metisazone (thiosemicarbazone) Anti-inflammatory sesterterpene | 1910-68-5 
Synthetic (heterocyclic 
34 Moroxydine biguanidine) Influenza treatment 3731-59-7 
Synthetic (Derivative of A nonnucleoside reverse 178040- 
35 Opaviraline N-Aryl-a-amino acids) transcriptase inhibitor 94-3 
Synthetic (ethyl ester A highly selective inhibitor of 
prodrug of oseltamivir influenza virus encoded 196618- 
36 Oseltamivir carboxylate) neuraminidase 13-0 
Synthetic (novel Specific inhibitor of human 153168- 
37 Pleconaril compound) picornaviruses 05-9 
Synthetic (nucleotide Antiviral treatment, necessary in 1809249- 
38 Remdesivir analogue) Covid research 37-3 
Synthetic (macrocyclic 923604- 
39 Simeprevir compound) Treatment of HCV 59-5 
Synthetic (amantadine Antiviral medicine used to treat 53783-83- 
40 Tromantadine derivative) herpes simplex virus. 8 


Binding site identification 

The binding pocket residues were selected on the basis of the literature survey and the PDBsum database. The amino 
acid residues were Thr 190, Glu 166, His 163, Phe 140, Gly 143, Cys 145, Gln 189, His 41, His 172, Ser 144, Asn 
142, and Leu 141. 


Molecular docking 

The molecular interaction analysis of the SARS-CoV-2 main protease (6LU7) and the selected 40 antiviral drugs was 
carried out in multiple docking software; MOE and AutoDock Vina. 

The initial level of molecular docking was performed using MOE software (molecular operating environment, version 
2015.10. The target structure purification, preparation, protonation, and defining of the binding site were performed 
as pre-docking steps. The prepared target and ligand structures were subjected to molecular docking using default 
parameters [7], 

The second level of molecular docking was performed using Autodock vina software (v.1.2.0.), comparing the 
interaction affinity of drug molecules to the target. Swiss PDB viewer was used for target preparation and docking was 
carried out in Autodock Vina. 

Interaction analysis of the receptor-ligand complex was performed using discovery studio bio via 2017. Poses with the 
best binding affinity and H-bond interaction were selected for comparison. The interaction and binding affinity of 
ligands in MOE and Autodock Vina were compared. Ligands with the least binding affinity and the most favorable 
hydrogen bonds were selected and subjected to an interaction validation study molecular dynamics simulation [8]. 


Molecular dynamic simulation 

The MD simulation analysis was performed using GROMACS version 2020.1 MD package in Ubuntu 20.04.2 (AMD 
Ryzen 9 3900 x 12-core processor x 24). The topology file of protein and ligands were generated separately using the 
Charmm 36 force field. Solvation, lonization, Energy minimization, and Equilibration using NVT (300 K), NPT (1 atm) 
ensemble were performed prior to production MD run for 100 ns. 


Trajectory analysis 

The trajectory analysis of the molecular dynamics simulation result was performed in the GROMACS version 2020.1 
package. The RMSD, RMSF, SASA, Rg, and Hbond plots were computed using gmx rms, gmx rmsf, gmx sasa, gmx 
gyrate, and gmx hbond tools, respectively. VMD, Pymol, and Discovery Studio Visualizer were used for the result 
visualization and analysis [9], 


MMPBSA 

Molecular Mechanics Poisson-Boltzmann Surface Area Method (MM-PBSA) was used as a scoring function for 
estimating the free energies between biomolecular interactions. The binding free energy of a complex in 100 ns MD 
trajectory was calculated using the command-gmx mmpbsa. 

To compute the individual energy contributions of each residue to the MGAM-TW/WA interaction, the binding energy 
was further decomposed on a per-residue basis. The binding free energy of a protein-ligand complex in a solvent can 
be estimated by 

AGpbind=AGcomplex—[AGproteintAGiig], 
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Where AGcompiex is the total free energy of the complex, AGprotein, and AGiig denote the energies of isolate protein and 

ligand respectively. Conceptually, the MM-PBSA approach can be described as 
AGbind=AEgast AGsoi=AEvawt AEele+ AG polartAGnonpolar 

where AEgas is the gas free energy, that is the average molecular mechanics potential energy in a vacuum which 

involves van der Waals (AEvaw ) and electrostatic (AEeie) interactions; AGsoi represents the contribution to the 

solvation-free energy which includes polar solvation (AGpolar) and nonpolar solvation (AGnonpolar) energies [1°. 


RESULTS AND DISCUSSION 


The in silico analysis of the FDA approved antiviral drug molecules against the SARS-CoV-2 main protease (PDB ID: 
6LU7) was performed to implement the drug repurposing approach. To prevent the pandemic, marine based and 
synthetic antiviral drugs were repurposed against SARS CoV-2. This study was carried out to explore the efficacy of the 
considered drug molecules against the virus. The computational study comprises a chain of steps including data 
collection, interaction analysis through molecular docking, and the validation of molecular interaction stability through 
molecular dynamics simulation [41-45], 

Interaction analysis of the docked complexes obtained from MOE and Autodock vina software was carried out using 
Discovery studio visualizer. The ligand poses were ranked based on their binding affinity and the number of H-bond 
interactions (Table 4). 


Table 4. Dock result of the marine drug molecules against 6LU7 ranked by affinity. 


SI | Source Ligand Name H-bond interactions No. of H Affinity 
No bond (kcal/mol) Affinity 
interactions | Autodock (kcal/mol) 
Vina MOE 
; GLN 189, ASN 142, SER 144, LEU 141, 
1 PseudopterosinF HIS 163, CYS 145, GLU 166 7 -7.9 -4.31878 
. HIS 163, CYS 145, GLY 143, LEU 141, 
2 PseudopterosinG SER 144, ASN 142 6 -7.9 -4.82708 
3 PseudopterosinE GLN 189, LEU 141, HIS 163 3 -7.8 -3.74498 
: GLU 166, HIS 163, CYS 145, GLY 143, 
4 PseudopterosinA LEU 141, SER 144, ASN 142 7 -7.4 -4.7229 
5 PseudopterosinC fice Ore cael pe oe eee 5 -7.4 -4.24546 
6 Manoalide SER 144, GLY 143, CYS 145, GLU 166 4 -7.4 -5.01197 
: ASN 142, LEU 141, GLY 143, CYS 145, 
7 PseudopterosinD HIS 41, GLU 166 6 -7.3 -5.28855 
8 Pseudopterosin X GLY 143, SER 144, CYS 145, THR 26 4 -7.3 -4.84206 
. HIS 163, SER 144, ASN 142, GLY 143, 
9 PseudopterosinB GLU 166 5 -7 -4.38097 
: Azidothymidine GLY 143, CYS 145, SER 144, LEU 141, 
10 | Marine (zidovudine) HIS 163 5 -7 -3.61424 
11 PseudopterosinW GIN te ee He AS ER Ey HEED: 5 -6.9 -5.64598 
GLU 166 
12 Polyacetylenetriol | LEU 141, HIS 163, THR 190, GLU 166 4 -6.9 -4.33643 
13 Sarcophine CYS 145, HIS 41 2 -6.7 -3.42628 
14 Sarcophytolide SER 144, CYS 145, GLY 143 3 -6.1 -2.41924 
. : GLU 166, GLY 143, SER 144, CYS 145, 
15 Vidarabine ASN 142 5 -5.9 -3.7503 
: GLU 166, HIS 163, ASN 142, CYS 145, 
16 Cytarabine SER 144 5 -5.8 -3.49828 
17 Chloroquine HIS 164 1 -5.8 -3.20447 
18 Sarcophytol GLU 166 1 -5.8 -2.57192 
; CYS 145, SER 144, LEU 142, PHE 140, 
19 Acyclovir GLU 166, ARG 188 6 -5.7 -2.84831 
20 Avarol HIS 164 1 -8 -3.52298 
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GLU 166, SER 144, CYS 145, LEU 141, 


21 Remdesivir HIS 163, GLN 189 6 -8 -5.64723 

22 Simeprevir LEU 141, SER 144, CYS 145, GLU 166 4 -8 -5.47218 

23 Glecaprevir eee nm oie i ood et ) -7.7 -6.67256 

24 L-756423 GLN 189, CYS 145, HIS 41, GLU 166 4 -7.5 -6.08067 

25 Indinavir ASN 142, CYS 145, SER 144, GLY 143 4 -7.1 -5.80856 

26 Lopinavir GLN 189, CYS ys 143, HIS 41, GLU 5 74 4.614 

27 JE-2147 THR 26, CYS 145, GLU 166 3 -7 -5.56536 

28 Amprenavir THR 24, HIS 41, CYS 145, HIS 5 6.9 5.90962 

29 Pleconaril aE ee ir sa dail 5 6.7 5.67369 
. GLN 189, GLY 143, CYS 145, SER 144, 

30 Synthetic Entecavir LEU 141, PHE 140, GLU 166 7 -6.5 -2.62398 

31 Tomantading: [oe ee eel ase, 5 6.3 -2.48673 

GLU 166 
32 Hydroxychloroquine | ASN 142, LEU 141, CYS 145, GLN 189 4 -6.2 -3.92233 
: ASN 142, GLY 143, CYS 145, THR 26, GLU 

33 Adefovir 166, PHE 140 6 -6.1 -3.53678 

34 Anagrelide ASN 142 1 -6.1 -3.72704 

35 Opaviraline HIS 41, CYS 145, GLU 166 3 -6.1 -5.17043 

36 Oseltamivir CYS 145, ASN 142 2 -6.1 -3.39892 

37 Metisazone LEU 141,ASN 142, SER 144, CYS 145 4 -5.7 -2.62684 

38 Moroxydine LEU 141, ASN 142, GLU 166, PHE 140 4 -5.2 -2.84513 
eatore rs HIS 164, SER 144, HIS 163, PHE 140, 

39 Favipiravir GLU 166, ASN 142 6 -4.9 -3.27563 
HIS 163, LEU 141, ASN 142, PHE 140, 

40 Foscarnet GLU 166 5 -4.1 -2.69203 


The drug molecules from marine based and synthetic were screened and listed separately. All the exhibited 
interactions are with the binding site residues. Among the 40 considered drug molecules, pseudopterosins (F, G, E, A, 
and C) from marine sources and Remdesivir, Simeprevir, Glecaprevir, L-756423, and Indinavir from synthetic drugs 
show better H bond interactions and binding scores in both MOE and Autodock Vina software [46], 

Five top ranked drug molecules from each category, synthetic and marine sources were subjected to molecular 
dynamic simulation in GROMACS software to validate the interaction stability (Table 5). 

The MD simulation result shows both stable and unstable interactions. Few ligands exhibit interaction in the same 
amino acid residues as in the molecular docking result and few shows different interaction due to the energy variation 
and conformational geometry changes [47], 


Table 5. Molecular dynamic simulation result of selected 10 receptor-ligand complexes. 
H-bond interactions in No of H bond 


S| No} Ligand name : ; H-bond interaction after MD result 
docked complex interaction 
GLN 189, ASN 142, SER 
1 Pseudopterosin F | 144, LEU 141, HIS 163, f GLU 166 


CYS 145, GLU 166 
HIS 163, CYS 145, GLY 


2 PseudopterosinG |143, LEU 141, SER 144, 6 GLY 179, PHE 181, PHE 185 
ASN 142 
GLU 166, HIS 163, CYS 
3 Pseudopterosin A |145, GLY 143, LEU 141, 7 GLU 166, GLN 192, THR 190 


SER 144, ASN 142 
HIS 41, CYS 145, GLY 
143, SER 144, ASN 142 


4 Pseudopterosin C 
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5 Pseudopterosin E oil marr ee 3 THR 45, CYS 44 
GLU 166, SER 144, CYS 
6 Remdesivir 145, LEU 141, HIS 163, 6 - 
GLN 189 
: GLU 166, ASN 142, CYS 
it Glecaprevir 145, GLY 143, HIS 44 5 GLN 189, GLN 192, THR 190 
8 Sianrevik LEU 141, SER 144, CYS A CYS 145, GLU 166, SER 144, GLY 143, HIS 
P 145, GLU 166 41 
9 Indinavir GLU 166 ASN 142 2 GLU 166 GLN 189 
10 L-752463 MET 49 1 - 


After molecular dynamics simulation, the trajectories were analyzed to understand the spatial fluctuations of protein. 
Among the subjected molecules, only pseudopterosin F, pseudopterosin G, pseudopterosin A, pseudopterosin E, 
Glecaprevir, Simeprevir, and Indinavir exhibit interaction after MD simulation. The drug molecules pseudopterosin F, 
pseudopterosin A, Simeprevir, and Indinavir show stable interactions in the GLU 166 residue. Molecular interaction is 
absent in pseudopterosin C and Remdesivir complexes. All the receptor-ligand complexes considered for the MDS were 
subjected to MMPBSA calculation and the result was analyzed (Table 6). The compound’s favorable binding energies 
may indicate their potency as inhibitors. The least binding affinity constitutes the strongest binding [48]. 

The binding energy result of the protein-ligand complexes ranges from 10.54 to -6.04 kcal/mol. Among the drug 
molecules, pseudopterosin A and Simeprevir show better binding energy -3.69 and -6.04 (kcal/mol) respectively. 
pseudopterosin A is a marine diterpene glycoside found in Caribbean Sea whips, pseudopterogoria elisabethae. The 
pseudopterosin group consists of thirty one structurally unique derivatives, that possess various functional properties 
such as antiviral, anti-inflammatory, antibacterial, anticancer, wound-healing, analgesic. Pseudopterosin A is one of 
the most potential and extensively studied pseudopterosin as they exhibit various mechanisms of action including the 
stabilization of cell membrane, inhibition of phagocytosis, intracellular calcium alteration, and act as a 
neuromodulatory agent. 


Table 6. The MMPBSA result of the selected drug compounds with their 2D structure. 


S| No Source Ligand name Binding energy Structure 
(kcal/mol) 
1 Marine Pseudopterosin F 10.54 
On 


2 Pseudopterosin G -0.71 
e- 
oe 
af 
I 
3 Pseudopterosin A -3.69 Ue 
be 
wane 
iaeer 
4 Pseudopterosin E 8.36 
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5 Synthetic Glecaprevir -1.80 
6 Simeprevir -6.04 
7 Indinavir 0.29 


Pseudopterosin A exhibits h-bond interactions in GLU 166, HIS 163, CYS 145, GLY 143, LEU 141, SER 144, and ASN 
142 during molecular interaction analysis from which GLU 166, GLN 192, and THR 190 residues were found to be 
stable after the MD simulation. Here the GLU 166 residue remains stable even after the 100 ns MD simulation (Figure 
1). In the MD trajectory analysis of the 6LU7-psuedoteridopsin A complex, the mean RMSD values in the backbone 
were 0.574 nm, and the RMSD values varied between 0.9788 nm and 0.0004 nm. The graph shows an observable 
deviation up to 63 ns and an equilibrated state from 63 ns to 95 ns with an average of 0.870 nm. The RMSF implies 
that the regions from residue index 274-279 show fluctuation between 0.82 nm to 1.085 nm, residue index 195-197 
show fluctuation between 0.796 nm to 0.925 nm, and residue index 72-78 show fluctuation between 0.799 nm to 
0.817 nm. The RoG of the complex fluctuates between 2.37 nm to 2.15 nm with an average of 2.23 nm. The 6LU7- 
psuedoteridopsin A complex has an average SASA value of 151.2239 nm? and ranges from 161.392 nm? to 141.617 


nm2 (Figure 2). 


Figure 1. Molecular interaction between psuedoteridopsin A and SARS-CoV-2 main protease (6LU7). (A) Interaction 
of psuedoteridopsin A-6LU7 complex after molecular docking; (B) Binding pose of docked complex, 
psuedoteridopsin A-6LU7 in hydrophobic surface view; (C) H-bond residues of docked complex, psuedoteridopsin A- 
6LU7; (D) Interaction of psuedoteridopsin A-6LU7 complex after MD simulation; (E) MD result of psuedoteridopsin A- 
6LU7 complex in hydrophobic surface view; (F) H-bond residues of psuedoteridopsin A-6LU7 complex after MD 


simulation. 
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Figure 2. Trajectory plots of psuedoteridopsin A-6LU7 MD simulation. (A) RMSD plot of psuedoteridopsin A-6LU7 MD 
simulation; (B) RMSF plot psuedoteridopsin A-6LU7 MD simulation; (C) Radius of gyration plot of psuedoteridopsin A- 
6LU7 MD simulation; (D) H-bond distribution plot of psuedoteridopsin A-6LU7 MD simulation. 
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Simeprevir is an oral drug that acts directly to inhibit hepatitis C virus protease (HCV NS3/4A) and is used in HCV 
infection in adults. As they are macrocyclic compounds, the structure itself helps to enhance the affinity and 
selectivity for speedy association and slow dissociation to the receptor molecule via non-covalent binding. This 
antiviral drug is highly efficacious and safe, with fewer side effects. The molecular interaction analysis exhibits h-bond 
interactions in LEU 141, SER 144, CYS 145, GLU 166, and molecular dynamic simulation exhibits h-bond interactions 
in CYS 145, GLU 166, SER 144, GLY 143, and HIS 41 residues. The residues CYS 145, GLU 166, and SER 144 present 
in molecular docking remain stable after MD simulation (Figure 3). Here most of the residues involve in receptor- 


ligand binding are stable [19]. 


Figure 3. Molecular interaction between Simeprevir and SARS-CoV-2 main protease (6LU7). (A) Interaction of 
Simeprevir-GLU7 complex after molecular docking; (B) Binding pose of docked complex, Simeprevir-6LU7 in 
hydrophobic surface view; (C) H-bond residues of docked complex, Simeprevir-6LU7; (D) Interaction of Simeprevir- 
6LU7 complex after MD simulation; (E) MD result of Simeprevir-6LU7 complex in hydrophobic surface view. (F) H-bond 
residues of Simeprevir-6LU7 complex after MD simulation. 


In MD trajectory analysis, the mean RMSD values in the backbone of the 6LU7-Simeprevir complex were 0.211 nm and 
the RMSD values varied between 0.3675 nm and 0.0005 nm. After the initial deviation, the RMSD value attains 
equilibrium from 50 ns to 75 ns with an average of 0.180 nm which is so close to the total average, followed by a 
slight slope in the plot to attain the maximum value of 0.3675 nm [29], 

The RMSF of the complex implies that both terminal ends showed higher fluctuations as they are highly dynamic in 
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nature. The highly fluctuating region is from residue index 46 to 56 with RMSF values ranging from 0194 nm to 0.3 
nm. Other residues showing high RMSF values are residues 169, 168, and 191 with values of 0.278 nm, 0.238 nm, 
and 0.225 nm respectively. The RoG fluctuates between 2.31 nm to 2.19 nm with an average of 2.25 nm. For the 
6LU7-psuedoteridopsin The SASA value for the complex has an average value of 152.8134 nm2 and ranges from 
161.319 nm? to 144.579 nm? (Figure 4). 


Figure 4. Trajectory plots of Simeprevir-6LU7 MD simulation. (A) RMSD plot of Simeprevir-6LU7 MD simulation; (B) 
RMSF plot Simeprevir-6LU7 MD simulation; (C) Radius of Gyration plot of Simeprevir-6LU7 MD simulation; (D) H-bond 
distribution plot of Simeprevir-GLU7 MD simulation. 
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CONCLUSION 


The current study is an in-silico approach to repurpose the commercially antiviral drugs against SARS-CoV-2. 
Molecular interaction analysis and further MD validation of a set of 40 drug molecules and the SARS-CoV-2 main 
protease were performed. The Binding free energy of the top-ranked receptor-ligand complexes was calculated using 
MMPBSA. In the study, pseudopterosins from marine sources (Caribbean Sea whip-pseudopterogoria elisabethae) 
shows noticeable interaction with the target molecule. The molecular dynamic simulation and MMPBSA results imply 
that pseudopterosin A and simeprevir show promising interactions and binding affinity against the receptor. The 
potentiality of these drug molecules against the virus brings new hope for COVID-19 medications. Further clinical 
analysis is needed to validate the activity. 
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